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I 1 Abstract 

o 

5_( ■ The cosmological expansion of the universe affects the behaviour of 

iy>j ' all physical systems and, in the case of gravitationally bound ones, could 

correspond to or mimic a time dependent Newton's constant. Here we 
discuss the case of a locally spherical mass distribution embedded in a 
generic Robertson Walker universe. Choosing the most appropriate metric 
tensor for the problem and assuming that the local time scale is much 
QO ' much lower than the cosmic one, we show that G is practically unaffected 

QO , thus leaving the absolute magnitude of type la supernovae unaltered at 

QO ■ all epochs. 

q ' 1 Introduction 

00 ' 

O 

An open and since long time debated problem is the one of the effect of the 
cosmic expansion on the behavior of gravitating systems on a more or less local 
^ ' scale (a few examples out of the wide litterature are [1] [2] [3] [4] [5] [6]). Of 

the various ways this problem can affect the understanding of our universe, we 
would like to focus on the collapse of a white dwarf leading to a supernova of 
type la [7] [8] (Snla). Snla's are particularly important because they have led 
to the discovery of the accelerated expansion of the universe. The fundamental 
feature of Snla's for the monitoring of the expansion is their being considered 
as "standard candles". The latter conviction is based on the fact that the 
peak luminosity of this kind of supernova is proportional to the mass of nickel 
synthesized during the collapse, which, to a good approximation, is a fixed 
fraction of the Chandrasekhar mass [9]. The Chandrasekhar mass, in its turn, 
is proportional to G~ 3 / 2 [10] where G is the universal coupling constant between 
matter and geometry (Newton's constant). 
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There are alternative theories to General Relativity (GR) for which G is no 
constant at all, but rather it depends on cosmic time (examples are [11] [12] [13] 
[14] [15] [16] [17]). In that case, provided one has the explicit time dependence 
of " G" , the absolute luminosity to be attributed to an Snla should appropriately 
be corrected on the basis of the corresponding redshift z (distance in time) in 
order to separate the contribution of the expansion from the one of the changing 
G. 

Even in the case of GR (possibly "extended" [18]), however, one could sur- 
mise an influence of the expansion of the universe appearing in the form of 
an effective G e ff, i.e. something "dressing" the universal G with a time- or 
z-dcpcndcncc of the local gravitational interaction. 

In any case the presence of a time (z) dependent G e // (whatever be the 
origin of the dependence) would imply a correction on the absolute magnitude 
M (z) of an Snla, given by the formula [19]: 

M(0)-M(,) = ^log^klM 

4 G eff (0) 

Considering the relevance of this problem, in the present letter we shall fo- 
cus on the supernova luminosity in any Robertson Walker (RW) spacely flat 
universe, i.e. any expanding universe in which the gravitational interaction is 
expressed by the geometry of a four dimensional manifold with a global symme- 
try reducing the number of independent functions, in the average cosmic metric 
tensor, to the only scale factor a. 

We shall start from an ansatz for the metric describing the situation around 
a local spherical distribution of mass in an expanding RW universe. Then 
exploiting the fact that a supernova explosion is an extremely short event on 
the scale of cosmic times we shall verify that the implied G value is in any case 
the universal one. 



2 A spherically symmetric perturbation in a Robert- 
son Walker universe 

Let us start from the remark that at the scale of a stellar system the effect of the 
cosmic curvature appears only as a very tiny perturbation of the usual stationary 
state solutions of general relativity. Considering for instance the example of a 
spherically distributed bunch of matter in a background RW universe with zero 
space curvature, we may guess and assume that locally the induced metric is 
essentially 

ds 2 = (1 - / (r, r)) dr 2 - a 2 (r) ( - d f. - + r 2 d9 2 + r 2 sin 2 6d^ 2 ) (1) 

V 1 - h (r, t) J 

where a is the cosmological scale factor (adimensional) and an extremely weak 
dependence of / and h on time is expected. Our r is measured in meters. 
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Far away from the local source it must be /, ft, — > so that the pure RW 
metric is recovered. On the other side if it were a = a = constant we could 
incorporate its value in the definition of r and from the field equations we would 
obtain the Schwarzschild solution f = h = 2Gm/c 2 r = 2/i/r. 

A metric like the one in (f) belongs to a general class of exact solutions of the 
Einstein equations discussed for instance in [20] . When applyed to a cosmologi- 
cal problem it has a partial correspondence with the early work of McVittie [1] , 
many times reconsidered and discussed afterwards (see for instance [21] [22]). 
Actually, however, he analyzed a spherical massive source immersed in a homo- 
geneous fluid with which it interacted gravitationally. Our case (and specific 
treatment) is different because for us the cosmic fluid is implicitly unaffected by 
the local spherical source, which assumtion is indeed realistic being the cosmic 
fluid itself the product of an average over the whole universe. What we want to 
elucidate is the inverse, i.e. the effect of the global (average) curvature on the 
local system. 

Without loosing generality we may restyle (1) evidencing an arbitrary mo- 
ment in cosmic time. The corresponding a (to) — an (the value of the scale 
parameter at the chosen moment) may be absorbed into a rescaling of r, and 
we introduce the new renormalized scale factor a (r, t ) = a (r) /a (tq). 

From the line element written in this form we can compute the Einstein 
tensor G M „, assuming that f,h«l which is an obvious situation when black 
holes are not implied (actually, in "ordinary" conditions, /, h <~ 10~ 6 ): 

Goo = 3— lh'+-)--h 

a. 1 ra z \ r J a 

G rr = X - + +aaf+(2aa + a 2 )(f + h-l) 

Gg = (2aa + a 2 ) r 2 (f - 1) + X -o?r 2 h (2) 

+aar 2 (j + ^hj + \r 2 f" + \r (/' + b!) 

G u = G ee sin 2 6 
r a h 

^Or — —J 

a r 

Primes mean partial derivatives with respect to r and dots correspond to partial 
t derivatives. 

One further assumption is that the source term in the Einstein equations 
can be decomposed into a cosmic contribution T M „ (peculiar to the given model 
one is considering) and a local contribution 1^, so that the equations can be 
written: 

The two contributions are necessarily very different from one another and the 
second is locally much bigger than the first. Having assumed, on the average, 
a RW symmetry, we can verify, by direct inspection of (2), that the terms 
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appearing in the expression of G^y not containing / and h correspond to the 
cosmic source, so that the remaining ones pertain to the local source. If we 
consider the situation outside the local distribution of matter we are in " vacuo" , 
so that the additional terms of (2) give rise to the equations 



J_ h 



ti + - + -h = 
r I a 



\ (/' + 7) + adi f + (2«a + « 2 ) (/ + h) = 

(2aa + a 2 ) r 2 f + ]-a 2 r 2 h + aar 2 (f + H) (3) 



+ \r 2 f" + \r{f' + h')=Q 
a r 

The assumed weak dependence on time allows for a low order power series 
development around the local (in time) expression of the functions. Introducing 
t = t — to as the time variable we write: 

a(t) = l + H Q t-^q H 2 t 2 + ... 

/(r,t) = f (r) + h (r) t+/ 2 (r) t 2 + ... (4) 
h(r,t) = ho(r) + h 1 (r)t+h 2 (r)t 2 + ... 

Just to fix the orders of magnitude let us remark that: 

f,h < 1CT 6 at most 
H ~ HT 18 s- 1 =~ 10" 26 m- 1 

The Hubble parameter correction can produce a contribution comparable to the 
one from / and h in times of the order of ~ 10 5 years. 

Introducing (4) into (3) and fixing the attention on the zeroth order in t we 
obtain: 

^ + ^ + HQhl=0 

7 + ^ + H o 0- ~ 2 *>) (/o + *»<>) + ff o/i = 
/i 2 r 2 + i/ r 2 (/1 + ^1) + y /o + «o (1 - 2 9o) r 2 /o (5) 



+1 (K + /o) = 

Ho/5-- = 

r 



4 



The system (5) is easily solved, starting from the ansatz: 

fo = 2*. 
r 



The solution is 



fo = 


2» 

r 








h = 










fi = 


A -2(i 
H r 3 


H --iTo(l 

r 


- 2q ) i 


(2(1 + A tt2 \ 
{ r + ff oV) 


hi = 


-2H » 
r 








h 2 = 


3 (A 

-A-2-v 


)+7^(i- 


2q ) + 


^ 7 -H + rt iH*(l- 



(6) 



where A is an integration constant. 

In terms of the full functions we would write 



/«,t, : 2^-(^ + ^ + H (l-2 ft )f^ + ^r))t + .. 



h(r,t) = -+H^(ir-2H -t 

r r 



+ (4f (7 - m) + 7^o 2 (1 - 2«) + ^ + ^ (1 - 2ft)) t 2 + 

In a flat, static background it would be H — and the solution should coincide 
with Schwarzschild's, so that it must be 

A = 2(i 

and 



/ (r, t) = 2^ - (i ^ + (1 - 2^0) + flgr 



t+... 



/»(r,t) = ^+H V-2ffo-t 

r r 



+^(rtf 2 (l-2 go ) + ^-4^ t 2 + ... 



We see that the expansion of the universe shows up in a "decoupling" of h 
from / already at the zero order. However this correction is extremely small. 
In fact the zero order for h is: 

2 ^ fi , H o 



r \ 2 

and the difference from the typical Schwarzschild term appears (1% correction) 
at distances ~ 10 25 m = 10 3 Mpc! 
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3 Behavior of test particles 



In order to investigate further on the possible local effects of the expanding uni- 
verse it is convenient to write down the geodesies for the metric in (1). Directly 
exploiting the spherical symmetry, which is not spoiled by the expansion, we 
may fix 6 = tt/2, so that the three remaining independent equations for the 
motion of a test particle are 



£l + ad£ + ad£ r 2 ( ^) - f - ( —) \ f'- — 

ds 2 \ds J \ ds J 2 V ds J ds ds 




*o» + o« (* + /)! 1^1 -oaf* 



d 2 r ^adrdr 
ds 2 a ds ds 



d0 
ds 



2a 2 \ds J 



j^dr dr h' ( dr 
ds ds 2 \ds 

d 2 <j) 2 dr d<f> ^adipdr 
ds 2 r ds ds a ds ds 
Eq. (7) is easily integrated to yield: 



hr ~~~ 



L 



ds 



being L a constant of the motion (angular momentum) . 
Using (8) the remaining two equations read: 



+«*»■» (WHS* 
•|(£) 2 + /'£f -(la 2 + aa(h + f) 



ds 2 



r - + 2 



a dr dr 



t (dr\ — hfkLfkL — hL (dr\ — D 
2a 2 \ ds) IL ds ds 2 \ ds) U 

We may also directly obtain dr/ds from the line element (1): 



- r (1 + h) 



(l-/)(l-/i) fdr\ 2 1-h 



a* 



ds J 



(1-h) 



(7) 



(8) 



(9) 



(10) 



The last three equations are what is needed to pursue any further analysis 
of the behaviour of matter in our system. 
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3.1 A special case: a fixed object 



We consider now a little object maintained at a fixed position r = R in the 
field (by means of some appropriate force); the angular momentum will then be 
L = 0. In this case equations (9) become: 



\r-=R (dry 
2a 2 \dsJ 



ds 2 2 V ds ) VXJ /I 1 \ 



where g is the four-vector representing the force per unit mass acting upon our 
object. 

At this point let us introduce the approximate solutions for a, /, h and keep 
the order in t: 

R rr n n q 5 \ / dr 



±l+H J2 q R+^ 



dg2 , ^y-*** . 2 "U . - H 2 R + 2R ) [ ds \ -00 



fi ( dr\ 



The first equation tells us that, for r = R, jj^ = K (R) also is a constant (to 
say it better: it is a function of R). 

This result is then introduced into the second equation, thus yielding (now 
g (i?) = c 2 -^rp^y and the explicit expression for fi has been restored) : 

-^r=9{R) (12) 

where g is in practice the local gravitational acceleration and of course (12) is 
Newton's law. 

What is indeed remarkable is that the coupling constant is the universal G, 
unaffected by the expansion of the universe. A phenomenon which happens in 
a short time, such as the implosion of a star due to a gravitational collapse, will 
be determined and measured by the same value of G at all epochs. This is of 
paramount importance, as we noticed before, for the evaluation of the apparent 
luminosity of type la supernovae, which depends on G~ 3//4 . 



4 Conclusion 

We have studied the case of a spherical distribution of matter embedded in a 
RW spacely fiat universe. Assuming that the global symmetry and behaviour of 
the universe are not affected by the local source, we have analyzed the influence 
of the cosmic expansion on the behaviour of the local system. What we have 
seen is that, for sufficiently small time spans (zeroth order expansion in time 
of the functions in the metric) the three-dimensional force acting upon a static 
test mass is indeed expressed by Newton's law without any epoch dependent 
renormalization of G. In this way the peak luminosity of Snla's, which depends 
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on G~ 3 / 2 and is due to a very rapid phenomenon, stays the same at all cosmic 
times. This result turns out to be valid for any RW type universe, i.e. for 
all theories where gravity is described by the curvature of a four-dimensional 
Riemannian manifold with isotropic homogeneous (in the mean) sources and flat 
three-dimensional space submanifold. This is true for instance for the Friedmann 
Robertson Walker, the A cold dark matter, and the cosmic defect theory [18]. 
In fact the result we have obtained could be anticipated considering that the 
gravitational interaction in an Einstein-type theory is always mediated by the 
unique coupling constant between matter and geometry, no matters where you 
are in cosmic time. Any cosmological effect shows up only on space and time 
scales far bigger than the ones of a local system. 
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